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Reduced wavenumber models of turbulence, shell models, show cascade processes and anomalous 
scaling of correlators which might be analogous to what is observed in Navier-Stokes (N-S) turbu- 
lence. The scaling properties of the shell models depends on the specific symmetries and invariants 
of the models. A new shell model is investigated. It is argued that this model might have a closer 
resemblance than the standard GOY model to the N-S turbulence. The new shell model coincides 
with the Sabra model proposed by L'vov et al. for a specific choice of the free parameters of 
their model. The model has for this choice of parameters, besides the energy and the 'helicity', a 
. cubic inviscid invariant. 
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The connection between Kolmogorovs 1941 (K41) scaling theory of turbulence and the N-S equation is through 
the four-fifth law ||. The four- fifth law is one of the few rigorous result regarding 3D hydrodynamic turbulence, 
connecting the second and third order longitudinal velocity structure functions. The equation is not closed, but in 
the inertial range, rj <C r <C L, r\ being the inner (Kolmogorov) scale and L the outer (integral) scale, we have 
S^(r) = — (4/5)er, where e is the mean energy dissipation. This implies £(3) = 1, where £(p) is defined from the 
qq ' scaling of the longitudinal structure functions, S p (r) = {5v\\(r) p ) ~ r^ p >. The scaling of all other structure functions 
, are observed to deviate from the K41 prediction, which from dimensional counting would imply = p/3. The 
deviation from the K41 prediction is called anomalous scaling exponents, referring to the intermittent nature of the 
energy dissipation, where the energy is inhomogeneously dissipated. 

As was noted recently || the inviscid conservation of helicity leads to another exact scaling law for a third order 
correlator associated with the flux of helicity, provided an inertial range for helicity flux exists B . Similarly for any 
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inviscid invariant an exact scaling law can be derived for the correlators associated with the spectral flux of such a 
quantity provided there is an inertial range separating sources and sinks for this quantity. This is the case for the new 

O ■ shell model, where the existence of a third order inviscid invariant leads to a scaling law for a fourth order correlator. 
, This has been argued to be a calculation of an anomalous scaling exponent ||ic|| . This is, however, not an anomalous 
scaling exponent in the afore mentioned sense, it would still be present in the case of an (imagined) homogeneous 
non-intermittent flow. On the contrary, the scaling of this specific fourth order correlator will be 'normal' in the same 

C3 ' sense as £(3) = li from the four-fifth law. 

The resent interest in shell models is mainly that they show numerically the same type of intermittent behavior 
as seen in 3D turbulence [|| . The simplicity of the shell models makes it possible to calculate the anomalous scaling 
exponents with high accuracy, which is still an undo-able task for the N-S equation, so in this sense shell models 
might prove useful as the starting point for exact results regarding calculating scaling exponents. The introduction 
of the Sabra model Q as superior to the GOY model § was motivated by improvements in accuracy with respects 
to numerical determination of scaling exponents. Here the emphasis will be on the similarity with the N-S equation 
and the cascade properties of the third order invariant. In fact the model will not exhibit a cascade of the third order 
quantity, which is easily seen from a scaling argument and confirmed in a numerical simulation. Thus it will be shown 
that the scaling exponent obtained in J10| is irrelevant. 



II. THE GOY MODEL 



The GOY model has built in the K41 scaling in the sense that the K41 scaling is a fixed point of the model. 
Furthermore, it has an unfortunate modulus(3) symmetry in shell numbers which has no resemblance in the N-S 
equation and which makes a precise numerical determination of scaling exponents difficult [6 10 1. The modulus(3) 
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symmetry, furthermore , introduces artificial long range (in /c-space) correlations with no analogs in the N-S equation. 
The GOY model has 2 inviscid invariants, the energy and a second non-positive definite quantity dimensionally 
equivalent to the helicity in 3D N-S turbulence. This 'helicity' only vaguely resembles the helicity in the N-S fluid, 
and it has been argued that it leads to an anomalous scaling behavior of the GOY model different from the mechanisms 
for intermittency in N-S turbulence p]]. A review of the main differences between the GOY model and the Sabra 
model is given in Ref. JjJ. 

III. THE NEW MODEL 

The shell model, defined in the following, can be regarded as a special case of the Sabra model introduced by L'Vov 
et al. It has the same two quadratic inviscid invariants, energy and 'helicity' as the GOY model. Furthermore it 
has one cubic inviscid invariant. The energy is the only positive invariant. Like for the N-S equation, and in contrast 
to the GOY model, the K41 scaling is not a fixed point of the model. 

The model is in the usual way defined by a set of exponentially spaced one-dimensional wavenumbers k n = kgX 71 for 
which we have u n as the complex shell velocity for shell n, (n = 1, N), The form of the governing equation for the new 
model is motivated from two demands. Firstly, the momenta involved in the triad interactions must add up to zero 
like in the N-S equation. Secondly, the complex conjugations involved in the non-linear term should be the same as 
for the terms in the spectral N-S equation involving triades for which the moduli of the wave vectors fall within three 
consecutive shells. This together with the usual construction of local interactions in fc-space, inviscid conservation of 
energy and fulfillment of Liouvilles theorem gives the equation of motion for the shell velocities, 

(d/dt + vk^)u n = 

i(k n u* n+1 u n+2 - ek n - 1 u* n _ 1 u n+1 + (1 - e)k 

n — 2 "> n — 1 ** n — 2 ) + fn (1) 

where v is the viscosity and /„ is the external forcing. The forcing would as in the GOY model typically be taken to 
be active only for some small wave numbers, ex. /„ = f5 n ^. 

Boundary conditions can be specified in the usual way by the assignment U—i = uq = «tv+i = un+2 = 0. 

The first requirement is fulfilled if the wave numbers k n are defined as a Fibonacci sequence, k n = k n —i + fc„_2- 
The choice of a Fibonacci sequence for the momenta leads to a model with the shell spacing uniquely being the golden 
ratio g, since for any choice of fci, k 2 (ki < k 2 ) we have k n /k n -\ — > g for n — > oo. So this corresponds to the usual 
definitions of the shell wave numbers with the golden ratio as shell spacing for ki =l,k 2 = (VE + l)/2 = j@. The 
golden ratio g = (1 + \/5)/2 plays a key role in the symmetries of shell models. 

With this formulation the shell spacing is not a free parameter of the shell model. However, using the definition by 
L'vov et al of k n = g n being a 'quasi-momentum' we shall keep the shell spacing, A, as a free parameter, k n = A™. 

If we interpret the momenta, k n , as representative of the modulus of the wavevector, k„, in 2D or 3D, the triangle 
inequality implies k n + fc„+i > k n+ 2 so the Fibonacci sequence corresponds in this sense to moduli of three parallel 
wavevectors. Note that for a shell spacing A > g (as the usual choice A = 2) the triangle inequality is violated. This 
means that we cannot interpret the usual shell model interactions as representative interactions between waves within 
three consecutive shells, since no such triplets of wavenumbers constitutes triangles. 

In order to give meaning to the notion of closing the triades, we define negative momenta, £;_„ = — k n , and assign 
the velocity, it_„ = it*, to these momenta. (The model still only has 2N degrees of freedom, represented by the N 
complex velocities). Note that ([!]) is also fulfilled for the negative momenta, which is why the prefactor must be V. 
With this notation we can rewrite (Q) as 

(d/dt + vk^+)u n = ik„ I(l,m\n)uiu m + /„, (2) 
h<k m 

where the sum is over positive and negative momenta, and all the dimensionless interaction coefficients have the simple 
form, 1(1, m;n) = I(l,m;n)S kn +ki+k m ,o with 1(1, m;n) = 8 n -%i$n-i,m - (eA)<5„_i,i<5„+i, m + ((1 - e) / 'X)S n +i,i^n+2,m- 
From this formulation it is clear why the complex conjugations in this model are exactly as in ([!]), they arise from 
the closing of the triades in the same way as for the spectral N-S equation. As noted in Q this is the main difference 
from the GOY model. 
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IV. INVISCID INVARIANTS 



It can easily be shown that there are only second order (quadratic) invariants of the form, ^ £ n u n U— n . The inviscid 
conservation of these quadratic invariants is obtained from, 

\n\<N 

i £, n k n I(l,m;n)uiu m u- n = 

ki~\-k m =k 7l 

i £ ?" _2fc n-2(l - e£ - (1 - e)C 2 )w„-2Wn-lti-n 
3<|n|<JV 

= 0. 

so exactly as for the GOY model we obtain the equation, 

1 - ef - (1 - e)e = (4) 
with the two solutions £ = 1 and £ = l/(e — 1). The first corresponds to energy conservation, with 

E = ^ -En = ( 5 ) 
and the second to 'helicity' conservation (for e < 1), with 

H = Y. H " = \ £ (tAt)"^-"- (6) 

M<jv v ' 

With the definition of negative momenta the only possible pth order invariants with terms u^.-.Ui must have the 
associated momenta summing to zero. Thus for any invariant, 

^2^lu h ...u jp + ... + i\u h ...ui p +c.c, (7) 

i 

the corresponding momentum vectors must add up to zero, kj 1 + ... + kj — ki t + ... + ki p = 0. This can easily be seen 
by differentiating (fit) with respect to time using (||). The only possible third order term is 

G = G n = X! C l (Wn-l«nU-(n+l) + C.C.) (8) 

n n 

By taking the derivative with respect to time using (||) we obtain a set of equations similar to (Q): 

e£ + A = 
(1 - e)£ 2 — A 2 = 
(1 - e)£ + eX = (9) 

These equations are fulfilled if 1 — e — e 2 = with £ = — A/e. The two solutions are, e = — g and e = 1/g. 
From a similar analysis it can also be seen that there are no invariants with p > 3. 



V. PARAMETER SPACE 



With e and A > 1 being the two free parameters of the model the parameter space is represented in figure 1 . The 
hyperbola is the curve A = 1/ 1 e — 1| corresponding to the dimensionally correct helicity, H — ^(— l) n fc n |w„| 2 . The 
dashed vertical lines corresponds to values where the third order quantity G is an inviscid invariant. The point, 
(e, A) = (l/g,g 2 ), marked with a square, has both the usual helicity and G = l)™5 3 "' ti n-i u n M -(n+i) conserved. 
The point (e, A) = (2 — g,g), marked with a diamond, does not have G conserved while the point (e, A) = (1/5,5), 
marked with a cross, has G conserved with £ = — g 2 and a helicity of the form, H = l)"/c I 1 / 2 |u„| 2 . The point 
(e, A) = (1/2,2), marked with a solid ball, is the point investigated by L'vov et al. corresponding to the values 
originally chosen for the GOY model. In the rest of this paper we will investigate the case (e, A) = (1/g, g 2 ) in order 
to discuss the role of the cubic invariant G. 
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VI. HAMILTONIAN STRUCTURE 



The Hamiltonian structure of the model as reported in ref. [jioj] can be observed, by change of variables, 



-n/2 



• i[AG/(-e)] 
Sv-„ 



(10) 



where G(v m , v— m ) is defined as a Hamiltonian with a density G n = (X\/ —e) n v n ^\v n v_^ n +i), which is local in wavenum- 
ber space. This relation more a curiosity than of practical importance. For shell models it would be more natural to 
have a Hamiltonian associated with the energy as in most conservative dynamical systems. The attempt to construct 
this has, however, until now not been fruitful |9j. 



VII. THE NON-LINEAR FLUXES 



The non-linear transfers of the invariants are defined as the currents, LT^ = (d/dt)^2 m<n E m , (v = f = 0), and 
correspondingly for H. They are, 



IE 



k n (D n+1 + (1 - e)D n /X) 



Yl% =k n (e-l)- n (D n+1 -D n /X), 

where D n = 2/m(u_(„_ 1 )M_„M n+ i). For e — 1/g the flux of G is, 

A? = (-AVeHAI^ - eDW/A - eDgO + d£>) 

where 

D ( r p = 2/m(u„ +2 w_(„ +1) 'u_ (rl _ 1 )M_ (n _2)) 
D { n 2) = 2Im(u n+2 u 2 _ n u_ (n ^ 1) ) 
D n ^ = 2Im(u n+2 u 2 _ (n+1) u n _i) 



(11) 
(12) 



(13) 
(14) 



VIII. STRUCTURE FUNCTIONS 



Scaling exponents are obtained from expressing the transfers in terms of structure functions. As found in Q 
the equation (^) is invariant under the rotation, u n — > exp(id n )u n where the phases are a Fibonacci sequence, 
Q n -2 + 9n-i = Qn- This symmetry is a trivial consequence of the construction of the model. The implication of the 



symmetry on the structure functions is that (uj 1 . 



"3pl 



exp[i(9 n + ... + ( 



jp 



)]<t 



"3pl 



Thus only structure functions 



with 9j 1 + ... + 8j p = can be non-zero. Since the phases fulfill the same relations as the associated momenta we can 
conclude that only structure functions where the associated momenta sum to zero are non-zero. The corresponding 
symmetry in the GOY model is Q n -i + #n-i + n = leading to the artificial slowly decaying modulus(3) correlation 
among distant shells. This is argued in |jj to make this model superior to the GOY model. 

The non-vanishing structure functions can easily be listed, thus we have the following 2nd, 3rd and 4th order 
structure functions, 



S. 



S 2 {n 
S 3 {n 

(0)/ 

4 (n,m 

Sf\n 
Sf\n 



= (u n U- n ) = (E n ) 

= 2Im(u n -lUnU_( n+ l)) = (D n ) 

= (KI 2 M 2 ) = (E n E m ) 

= 2Im(u n _ 2 Un-iu n+ iu_( n+2) ) = - (D^) 
= 2Im(u n - 2 u 2 _ n u n+ i) = - (D n 2) ) 
= 2Im(u n - 2 u 2 n _ 1 u_ {n+l) ) = -(D n 3) ) 



(15) 
(16) 

(17) 

(18) 

(19) 

(20) 
(21) 
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IX. EXACT SCALING RELATIONS 



The exact scaling relations corresponding to the four-fifth law simply states that the mean non-linear transfers 
jjE,h,G _ (j^e.HiG^ are independent of wave number within the respective inertial ranges in the high Reynolds 
number limit, 11^ = e, 11^ = 8 and 11^ = rj, where e,5,rj are the mean dissipations of E,H,G respectively. This is 



usually expressed in terms of structure functions. From (|ll|) and (16) we readily obtain, 



k n S 3 (n+l) + k n - 2 S 3 (n) = e (22) 
k n S 3 (n + 1) - k n ^S 3 (n) = (-iTk-^S. (23) 

in the inertial range, with the solution, 

fcn(l + A) 

The second term on the right hand side is an oscillatory term which is subleading in the scaling with k n in comparison 
to the first term. This term disappears in the case of a 'helicity free' forcing where 5 = and we recover the shell 
model correspondence to the 4/5th law expressed in terms of the third order structure function, S 3 (n) = e/k n . 
The equation similar to ( p2| ) for the transfer of the third order quantity, G, reads, 

kK-enXS^in + 1) + eSi X) (n)/A - eSf \n) + sf\n)\ = 

k 2 + a F i (n) = rj. (25) 

where rj is the mean dissipation of G, fifYi) denotes the square bracket on the left hand side, and a = log(— e) / log(A) = 
(wr+log(e))/ log(A). This was argued by L'vov et al. to establish a non-trivial calculation of a (sub-leading) scaling 

exponent, which in this notation reads, F 4 (n) ~ k~ 2 ~ a = k„ =>■ C(4) = 2+ (i7r + log(e))/ log(A). The imaginary part 
of the scaling exponent comes from the (—1)™ factor which can be trivially eliminated by reformulating the model in 
terms of new variables, v 2n = u 2n and v 2n +i = ~ u 2n +i^ as was done by L'vov et al. pC| . 



X. INERTIAL RANGES 



The validity of the exact scaling relations depends on the existence of inertial ranges separating the sources and 
sinks for the inviscid invariants solely associated with the non-linear fluxes of the invariants. The inertial range for the 
energy flux is determined by the Reynolds number Re, which we here just associate with the inverse of the viscosity 
v, assuming velocity at the integral scale to be of order unity. The Kolmogorov scale is in complete analogy to the 
K41 theory determined by balancing dissipation and non-linear flux, so we have Kg ~ (e/i/ 3 ) 1 / 4 growing as Re 3 / 4 . A 
similar analysis can be done for the analogous Kolmogorov scale for dissipation of helicity [W . Balancing dissipation 

and helicity input using u n - (e/fc„) 1/3 gives, K H - \f /(is 3 e 2 )] 1/T . Thus we get K H /K E - ^~ 3 / 7 + 3 / 4 = j, 9 / 28 -> 
for v — > 0. This means that for high Reynolds number flow the small scales will always exhibit non-helical flow. In the 
shell model the helicity changes sign due to the (—1)™ factor. The dissipation of positive helicity at even numbered 
shells and negative helicity at odd numbered shells will grow with wavenumber as - (-l) n fc 3 |w n | 2 - (-l) n k 7 n /3 , 
consequently the shell model will show strong odd-even oscillations of 11^ from balancing the positive and negative 
dissipations, the scaling will be determined by the dissipation |II^| ~ kj/ 3 for k n > K^. The situation for the cubic 
invariant G is different. We can again define a Kolmogorov scale for dissipation of G by equating the the dissipation 
and the input of G, 

vk 2 n G n ~ v{-l) n k 7 J 2 n[u n ^u n u*_ {n+l) ] ~ 77 (26) 

using u n ~ (e/fc„) 1/3 again gives, 

K G ~ [v/{ev)] 2 l\ (27) 

The ratio of dissipation scales is then Kg/ Ke ~ ^~ 2 / 5 + 3 / 4 = jy 7 / 20 , so as for the case of helicity the small scales will 
have no net non-linear flux of G. There is, however, a crucial difference between the dissipation of H and of G. The 
dissipation of helicity is forced to be of alternating signs whereas the dissipation of G can be of either sign at any shell. 
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Thus one should expect the mean dissipation of G n to vanish for k„ > Kg, and the dissipation would not determine 
the scaling of II„ . So how would the scaling of II„ be then? The non- linear flux 11^ is constituted of terms of the 
form k^J 2 Si{n) ~ fc„ ~ kj/ 6 . This leading scaling is killed by detailed cancellations between terms in an inertial 
range. One would then expect Kg to represent a de-correlation scale where the individual terms become independent 
and the scaling of |II^| becomes approximately that of the individual terms, which is kj/ 6 . This conjecture is tested 
in a numerical simulation. 



XI. NUMERICAL SIMULATION 



A simulation using a simple fourth order Runge-Kutta scheme has been performed with the parameter values, 
(e, A, v, N, f) = (1/g, g 2 , 10~ 9 , 20, (1 + i)5i^ n ). The anomalous scaling exponents are defined from ((ITf ) p / 3 ) ~ fen . 
The results for this model are show in figure 2 (diamonds) and coincides within the numerical accuracy of the simulation 
with the exponents found by L'Vov et al. for the Sabra model with (e, A) = (1/2, 2) (triangles). The non-linear fluxes 
for E, H, G are shown in figures 3, 4 and 5. The fluxes for H and G fluctuates between positive and negative values, 
so that only the absolute values are shown. The scaling indicated by the straight lines confirms the conjectures made 
in the previous section. However, even though the simulation is long enough and numerically accurate enough to 
determine the anomalous scaling exponents reliably, the non-linear flux II„ can only be determined using extreme 
numerical precision. If we conjecture that there is a constant flux of G of order unity (77 ~ 1) through the inertial 
range we must calculate 11^ as differences of fourth order correlators of order kl/ 6 which is about 10 6 at the end of 
the inertial range to obtain a constant of order unity. So in fact the graph in figure 5 does not reliably represent 11^, 
it is merely numerical noise. To see this consider the average |(II^)(£)| = | J Q U^(t)(1t\ as a function of t, sec figure 
6. The line is the curve a/ \fi which is expected for an independent random process. The standard deviation a of the 
process is of the order kj/ 6 . The main justification for studying shell models is the possibility of accurate numerical 
calculations of correlators and scaling exponents for high Reynolds number flow. Here we see that even the shell 
model can been pushed to the limit where the determination of correlators by numerical simulation is impractical. 



XII. SUMMARY 



To summarize, it has been argued why this model is a natural choice for a shell model of turbulence. For the 
choice of parameters which conserves both energy, dimensionally correct helicity and the third order quantity G the 
anomalous scaling exponents are the same as for the Sabra model which conserves energy and dimensionally correct 
helicity but not G. So even though the model has a Hamiltonian structure, the non-positive Hamiltonian, G, seems 
not relevant for determining the scaling properties of the model. It has been argued that there will not be an inertial 
range scaling regime for the fourth order correlator associated with the flux of G. However, due to the extreme 
numerical accuracy required to determine this correlator the numerical simulations presented here are not conclusive 
in the determination of this correlator. 
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FIGURE CAPTIONS 



Fig. 1 The parameter space (e, A) for the shell model. Energy is by construction always conserved. The hyperbola 
corresponds to conservation of dimensionally correct helicity H = ~}Z n {—^) n k n \u n \ 2 the horizontal line is the 
golden ratio shell spacing where the 'triangles' close, fc„_ 2 + fc n -i = k n . The vertical line represents parameters 
for which G = J2 n (— V e )™^[ M «-i M « u -(ri+i)] is an inviscid invariant. At the cross 'helicity' has the form 

# = E„(-i)"^ /2 KI 2 . 

Fig. 2 The anomalous scaling exponents ((p) defined from (n^/fc„) p / 3 <~ fcn^ p ' for (e, A) = (l/g,g 2 ) (diamonds) 
coincides within the numerical uncertainty with the values found by L'Vov at al. for (e, A) = (1/2, 2) (triangles). 

Fig. 3 The non-linear flux of energy as a function of k n . 

Fig. 4 The non-linear flux of helicity |II^ | as a function of k n . An inertial range with coexisting cascades of both energy 
and helicity is seen for the first few shells, after which the dissipation of helicity dominates the spectrum. The 
sign of the flux alternates for even and odd shells corresponding to dissipation of positive and negative helicity. 
The straight line indicates the scaling exponent 7/3 as is expected from the dissipation. 

Fig. 5 The non-linear flux of the third order quantity II^ as a function of k n . The numerics has not converged and 
the scaling exponents 7/6 indicated by the line comes from the standard deviation of the numerical value of the 
difference between correlators of the order kj/ 6 . 

Fig. 6 The average |(II^)(i)| = | J"' H^(t)(It\ as a function of t. The straight line has the slope —1/2 as for an 
independent random variable. This shows that the quantity shown in figure 5 is dominated by noise. Note that 
the sampling time is long enough to determine the anomalous scaling exponents. 
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